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Abstract 



T— I Poisson brackets provide tiie matiiematical structure required to identify the reversible contribution to dynamic phenomena in 
nonequilibrium thermodynamics. This mathematical structure is deeply linked to Lie groups and their Lie algebras. From the 
characterization of all the Lie groups associated with a given Lie algebra as quotients of a universal covering group, we obtain 
a natural classification of rheological models based on the concept of discrete reference states and, in particular, we find a clear- 
cut and deep distinction between viscoplasticity and viscoelasticity. The abstract ideas are illustrated by a naive toy model of 
crystal viscoplasticity, but similar kinetic models are also used for modeling the viscoplastic behavior of glasses. We discuss some 

0^ implications for coarse graining and statistical mechanics. 
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X. Which mathematics behind which physics? 

, The purpose of this work is to explore an elaborate math- 
ematical framework and its basic concepts in order to illumi- 
nate some physics and to develop a natural and general sce- 
nario for a range of phenomena. The underlying mathematics 
;is the abstract theory of Lie groups and, in particular, the corre- 
spondence between Lie groups and Lie algebras. Our concrete 
interest is in groups of space transformations. The physics asso- 
ciated with space transformations is the theory of deformation 
of materials, ranging from fluids to solids and encompassing 
all kinds of intermediate behavior. We are thus dealing with 
plasticity theory and rheology, viscoelasticity, viscoplasticity, 
and all that. We here make an attempt at classifying models 
^nd phenomena in rheology from a mathematical perspective. 
The key idea is that we identify discrete groups of equivalent 
reference states that evolve slowly in time, much slower than 
the structural variables describing deviations from the reference 
states. We lay the foundations for the kinetic theory of systems 
with reference states so that it becomes possible to bridge the 
gap between atomistic and phenomenological models of vis- 
coplasticity. 

[ Lie algebras are known to play an important role in the the- 
ory of complex fluids [19,1. Lie algebras allow us to introduce 
'Poisson brackets and Poisson brackets allow us to introduce re- 
versible dynamics (see Sec. |3]for details). With the help of a 
Poisson bracket, an energy function generates Hamiltonian dy- 
namics, which is the prototype of mechanistically controlled or 
reversible dynamics, to be separated from the irreversible "rest" 
in any thermodynamics. Nonequilibrium thermodynamics can 
take place only in structured spaces, more precisely, thermo- 
dynamic systems must possess a Poissonian structure. This 
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geometric structure is rooted in Lie algebras. More carefully, 
one should consider not only Lie algebras but also Lie groups, 
because one is usually working in infinite-dimensional spaces: 
Not every infinite-dimensional representation of a Lie algebra 
is guaranteed to exponentiate to a representation of some group 
(see §6.5 of Rossmann ll25ll '). The relevant Poisson brackets are 
actually obtained by a procedure known as Lie-Poisson reduc- 
tion from the canonical Poisson bracket on the cotangent bun- 
dle of a Lie group 1 15, 19], so that a group, or a representation 
of a group, is the true starting point for providing structure in 
nonequilibrium thermodynamics even though the final Poisson 
bracket is given in terms of the Lie algebra. 

Once one is aware of the fact that thermodynamics requires 
structured spaces associated with Lie groups and their Lie alge- 
bras, it is natural to look at all possible Lie groups associated 
with a given Lie algebra. These groups can actually be charac- 
terized very nicely ||8l|251]. There is a unique simply connected 
universal covering group Q, and all groups with the same Lie 
algebra are isomorphic to a quotient group of the form QIZ., 
where J3 is a finite or infinite discrete normal subgroup of Q 
(a subgroup J3 c ^ is normal if gzg^^ e Z for aW z & Z and 
g e 0). Actually, Z. must be contained in the center of the 
group 0, that is, all elements of Z commute with all elements 
of @ (or, gzg^^ - z; see §2.6 of Rossmann fz?]). What is the 
physical significance of the groups Z and QIZ in our applica- 
tion of group theory? Why can we interpret the elements of the 
discrete group Z as reference states? 

For illustration, let us consider the subgroup of unidirec- 
tional shear deformations of a simple cubic lattice in one of 
its principal directions, ^ = R, with addition as the binary 
group operation. For small shear deformations, we expect elas- 
tic behavior. For a unit shear deformation between neighbor- 
ing lattice planes, we reach a state that is fully equivalent to 
the undeformed state and, therefore, such a unit shear defor- 
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mation is plastic. The discrete normal subgroup of integers, 
J3 = 2, represents a set of equivalent reference states, and 
R/Z - Iq = [-1/2, l/2[ describes shear deformations with 
respect to a reference state. The topology of the Lie group 
R/Z (with the identification of -1/2 and 1/2) is that of a cir- 
cle, which is not simply connected. The occurrence of discrete 
subgroups is a general mathematical feature and it is natural to 
associate them with the equivalent reference configurations of 
a lattice. Only the simply connected universal covering group 
is free of reference states and hence provides the natural setting 
for viscoelastic liquids and other complex fluids. All smaller 
groups are associated with reference states and viscoplasticity. 

The detailed elaboration of the above picture is the purpose of 
the following sections, after compiling some basic knowledge 
from group theory (Sec.|2]) and nonequilibrium thermodynam- 
ics (Sec. [3]). By means of a toy example, we illustrate modeling 
on the full group (Sec. [Hi and on the quotient group (Sec. |5]) 
and observe that modeling on the quotient group is much sim- 
pler but requires some additional ingredients. A classification 
of rheological models resulting from the elaborated picture is 
offered in the final discussion (Sec.|6l). 

2. Some basic equations for groups 

The Lie group of interest in describing the response of mate- 
rials to stress is the group of smooth space transformations. If 
we consider a domain D c R^, the group of space transforma- 
tions is given by 



Q = {g : D ^ D\g smooth, g bijective). 



(1) 



with the composition of functions as the binary group operation. 
In the following sections, we mostly consider the subgroup of 
unidirectional shear transformations on D = R^, 
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where the binary operation of composition leads to the addition 
of the displacements A(y). The subgroup 0^ hence is commu- 
tative. The local shear strain is given by -/{y) = dA(y)/dy. For 
homogeneous shear deformations, we have A = yy with con- 
stant y. 

As small transformations can be visualized as effected by a 
velocity field acting over a short time, the Lie algebra associated 
with the group of space transformations on D provides the 
Poisson structure on the space of velocity fields on D lfl5i[l9ll . 
In the theory of complex fluids, the internal structure of the 
fluids is described by a vector space of functions defined on the 
same spatial domain as the velocity field. 



V = (v», 



D^R" 



(3) 



For the simplest case of Newtonian fluids, we have m = 2 for 
the number of structural variables because mass density and 
temperature are sufficient to describe the local state of the mov- 
ing fluid. For complex fluids, additional fields are needed to 
characterize the internal state of the fluid fl^]. To introduce 



structure on the space V one needs an action of Q on the dual 
of V, which is also a vector space. When Lie-Poisson reduc- 
tion is applied to the semidirect product of § and the dual of 
V, one obtains an extension of the Poisson structure from the 
space of velocity fields to the space of structural variables V 
ifislfigll . Typically, a space transformation g € acts on a real- 
or vector-valued function on D by transforming the arguments 
of the function. For real valued functions, the values at the 
transformed positions may be unchanged or multiplied by the 
Jacobian of the transformation, depending on whether we deal 
with a scalar or scalar density field. For vector or tensor func- 
tions, in addition, the components at the transformed positions 
may be mixed by proper tensor transformation laws. 

By introducing velocity fields and the structural variables (O 
in a region of space D, we have chosen to use an Eulerian de- 
scription of complex fluids. For a deeper understanding of the 
required Poisson brackets we can start with a Lagrangian de- 
scription of fluids, where the particle relabeling symmetry is 
the key to successful reduction to the smaller space of the Eu- 
lerian description (see ll26ll . Section 1.5 of 115 1, Appendix B.4 
of il9ll . and references therein). 

As mentioned before, the theory of complex fluids is based 
on functions defined on the same domain D as the group of 
space transformations These functions give the values of 
the structural variables at each point occupied by the fluid. An- 
other natural possibility to construct properly structured spaces 
employs functions defined on the group § itself, or some homo- 
morphic group. More precisely, we consider the vector space of 
integrable real-valued functions (the generalization to vector- 
valued functions is straightforward). 



/ = {/ : ^ ^ R). 



(4) 



The definition of / assumes a concept of integration and hence 
a measure on the Lie group ff. The natural choice is the Haar 
measure defined on all locally compact topological groups, 
which is invariant under the action of the group (see §5.2 of 
Rossmann ^23]). Up to a positive multiplicative constant, there 
exists a unique left-invariant and right-invariant Haar measure 
(in general, the two measures do not coincide). In particular, we 
are interested in the subset P c / of nonnegative functions with 
integral unity which can be interpreted as the set of probability 
densities for finding a particular space transformation or state 
of deformation of a (piece of) material. In kinetic theory, such 
probability densities are also referred to as configurational dis- 
tribution functions. In a stochastic interpretation, the evolution 
of the configurational distribution function can equivalently be 
considered as a stochastic process in the underlying Lie group 
Q (in particular, Fokker-Planck or diffusion equations for con- 
figurational distribution functions are equivalent to stochastic 
differential equations (H[3l). For local theories of the defor- 
mation behavior of materials, it may be sufficient to consider the 
subgroup of homogeneous space transformations (to be used 
separately at each point in space). 

Given a probability density p e P on ff, one can naturally 
define a probability density p^: on the quotient group Q/Z by 
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superimposing the contributions from all equivalence classes, 
Pzis)-YjP^^S)- (5) 

However, in general, there is no "periodic continuation" of a 
given probability density on because X may be an infinite 
group so that there arises a problem with normalization. More- 
over, one would need suitable "periodic boundary conditions" 
to obtain a smooth continuation. We denote the averages of a 
random variable f onQ or: QjX with respect to the probability 
density p or px as (/) or {f)z^ respectively. If a function / 
possesses the "periodicity property" f{zg) - fig) for all z e 2^ 
and g e Q, then we have (/) - {f)z- We here use the language 
of Fourier transforms because there exists a well-developed the- 
ory of harmonic analysis on locally compact groups, where the 
Fourier transform takes functions on a group to functions on the 
dual group, with particularly powerful results for commutative 
or compact groups. In the case of unidirectional shear transfor- 
mations of lattices, we deal with standard Fourier analysis and 
periodic solutions. 



3. Background from nonequilibrium thermodynamics 

Time-evolution equations for nonequilibrium systems have a 
well-defined structure in which reversible and irreversible con- 
tributions are identified separately. As pointed out before, the 
reversible contribution is generally assumed to be of the Hamil- 
tonian form and hence requires an underlying geometric struc- 
ture which reflects the idea that the reversible time evolution 
should be "under mechanistic control." The remaining irre- 
versible contribution is driven by the gradient of a nonequilib- 
rium entropy. We need a separate geometric structure to be ob- 
tained by entirely different arguments; however, for preserving 
symmetries under coarse graining, group representations play 
an important role also in the construction of the proper geomet- 
ric structure for generating the irreversible contribution to time 
evolution (see Section 6.1.6 of lfl9ll '). 

Our discussion is based on the GENERIC ("general equa- 
tion for the nonequilibrium reversible-irreversible coupling^ 
formulation of time-evolution for nonequilibrium systems IlZI 



dx 6E 6S 
— = L — + M — , 

dt 6x 6x 



(6) 



where x represents the set of independent variables required for 
a complete description of a given nonequilibrium system, E and 
S are the total energy and entropy expressed in terms of the vari- 
ables X, and L and M are certain linear operators, or matrices, 
which can also depend on x. The two contributions to the time 
evolution of x generated by the total energy E and the entropy 
S in Eq. (|6]) are the reversible and irreversible contributions, 
respectively. Because x typically contains position-dependent 
fields, such as the local mass, momentum and energy densi- 
ties of hydrodynamics, the state variables are usually labeled by 
continuous (position) labels in addition to discrete ones. A ma- 
trix multiplication, which can alternatively be considered as the 



application of a linear operator, hence implies not only summa- 
tions over discrete indices but also integrations over continuous 
labels, and 6/6x typically implies functional rather than partial 
derivatives. Equation (|6]l is supplemented by the complemen- 
tary degeneracy requirements 



6S 
L— =0 

6x 

and 

SE 
M — =0. 

5x 



(7) 



(8) 



The requirement that the entropy gradient 6S / 6x is in the null- 
space of L in Eq. O expresses the reversible nature of the L- 
contribution to the dynamics: the functional form of the entropy 
is such that it cannot be affected by the operator generating the 
reversible dynamics. The requirement that the energy gradi- 
ent 6E/6x is in the null-space of M in Eq. (O expresses the 
conservation of the total energy in a closed system by the M- 
contribution to the dynamics. 

Further general properties of L and M are discussed most 
conveniently in terms of the Poisson and dissipative brackets 



. 6A SB 
[A,B} = — L — , 
6x 6x 

6A 6B 
[A,B] = — M— , 

ox ox 



(9) 



(10) 



where A, B are sufliciently regular real-valued functions on the 
space of independent variables. In terms of these brackets, 
Eq. (|6]l and the chain rule lead to the following time-evolution 
equation of an arbitrary function A in terms of the two separate 
generators E and S , 



^ ={A,£) + [A,5]. 
dt 



(11) 



The further conditions for L can now be stated as the antisym- 
metry property 



{A,B) = -{B,A}, 

the product or Leibniz rule 

{AB,C}^A{B,C} + B{A,C], 

and the Jacobi identity 

{A, {B, C]] + {B, {C, A}) + {C, {A, B]] = 0, 



(12) 



(13) 



(14) 



where C is another arbitrary sufficiently regular real-valued 
function on the state space. The Jacobi identity (fT4l i. which is 
a highly restrictive condition for formulating proper reversible 
dynamics, expresses the invariance of Poisson brackets in the 
course of time (time-structure invariance). All these properties 
are well-known from the Poisson brackets of classical mechan- 
ics, and they capture the essence of reversible dynamics. 

Further properties of M can be formulated in terms of the 
symmetry condition 



[A,B] = [B,A], 



(15) 
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and the non-negativeness condition 
[A, A] >0. 



(16) 



This non-negativeness condition, together with the degeneracy 
requirement (|7]i, guarantees that the entropy is a nondecreasing 
function of time. 



dS 6S 6S 

= —M— = [5,5] >0. 
at 6x 5x 



(17) 



The properties (flST l and (fTST i imply the symmetry and the 
positive-semidefiniteness of M (for a more sophisticated dis- 
cussion of the Onsager-Casimir symmetry properties of M, see 
Sections 3.2.1 and 7.2.4 of Ottinger ['19']). From a physical 
point of view, M may be regarded as a friction matrix (actu- 
ally, often rates or inverse friction coefficients occur in M). 

4. Toy model on universal covering group 

We illustrate the general ideas of Lie algebras and quotient 
groups in the context of thermodynamic modeling for a toy 
model of crystal viscoplasticity OlOl] . We consider the unidi- 
rectional shear deformations (|2]i of a simple cubic crystal in 
one of the principal directions caused by shear stresses on the 
boundaries. The group @s of Eq. (|2]i is fully characterized by 
the unidirectional displacements A(3'). For homogeneous defor- 
mations, the group is isomorphic with the additive group of real 
numbers, y - [A(3')/3'] 6 R, and the set of equivalent refer- 
ence states is given by the discrete normal subgroup of integers 
J3 = 2. In this section, we model on the simply connected uni- 
versal covering group ^ = R of all shear deformations. This 
type of modeling on the full space is familiar from the theory of 
complex fluids, avoids the formulation of boundary conditions, 
and provides the most complete information about the evolution 
of reference states and relative deformations. Direct modeling 
on the quotient group R/Z is the topic of the subsequent sec- 
tion. 

Our toy model should not be confused with a realistic model 
of crystal viscoplasticity. It is well-known that the key to under- 
standing crystal deformation are dislocations. However, the dy- 
namics of dislocations involves the shearing of crystal domains 
and our toy model could hence offer a deeper understanding 

of the parameters (such as viscosity) in phenomenological the- 

I — I 

ories of dislocation dynamics [9]. Surprisingly, our model is 
also closely related to theories of shear transformation zones in 
metallic glass-forming liquids. For amorphous systems, peri- 
odic functions are used to model the average potential energy as 
an idealization of a much more complicated energy landscape 
with many minima on a wide range of energy levels 11211 . 



4.1. Thermodynamic formulation 

In order to formulate a toy model of crystal viscoplasticity, 
we introduce the momentum density field, which has a nonzero 
component only in the x-direction and depends only on the y- 
coordinate, g(y), the internal energy density e(y), and the proba- 
bility density p{y, y) over the shear deformations y between ad- 
jacent layers of a simple cubic lattice with spacing a. For every 



y, the configurational distribution function p{y, y) is a probabil- 
ity density. 

The total energy is given as the sum of kinetic, internal, and 
configurational energy contributions. 



E = Na- 



g{yY 



2p 



dy + 



p(y,y)^(y)dydy, (18) 



where is the number of atoms in a single displaced layer, Na^ 
is the area of this layer, and <1) is the potential energy associated 
with shearing a whole layer. For shear flows, a constant initial 
density p does not change with time and hence plays the role 
of a fixed parameter. The idea of equivalent reference states re- 
quires a periodic potential energy function O with period unity. 
In other words, the potential energy O is invariant under shear 
deformations y e Z. The prototypical example of such a poten- 
tial O is 



1 



2n^ (cos ny)^ + s ' 



with the force 



sin Iny 



2n [(cos Try)' + e]^ 



(19) 



(20) 



The parameter // has the dimensions of energy; Hy^ jl is the to- 
tal energy for small shear displacements y between two neigh- 
boring layers and small e. The parameter H hence has an ex- 
tensive character, H - HN, where H characterizes the energy 
per atom. The shear modulus is given by G = Hja^ so that, for 
known lattice spacing a, H can be obtained from a mechanical 
measurement in the elastic regime. The larger the number of 
atoms in the layer, the larger is the energy required to achieve 
shear deformations, and the smaller are the thermal shear fluc- 
tuations for a given value of k^T . The parameter e allows us to 
make the energy barrier between the different minima for each 
atom very high. In particular, the barrier becomes insurmount- 
able for e = 0. The control parameter e hence is the key to 
identifying well-defined, long-living reference states and defor- 
mations with respect to those. 

The total entropy is given as the sum of a thermal and a con- 
figurational contribution. 



S = Na- 



J s(e(y))dy - J J p(y,y)\np(y,y)dydy, (21) 



where s{e) is the density of the entropy associated with inter- 
nal energy and ke is Boltzmann's constant. The configurational 
entropy contribution has the well-known Boltzmann form. 

Before we can evaluate the functional derivatives of energy 
and entropy, we need to define the relevant linear spaces more 
carefully. In the spirit of Lie-Poisson reduction, we consider 
the functional derivatives of observables as elements of an un- 
derlying Lie algebra, and variations of the independent fields as 
elements of the dual Lie algebra. For every y, the normalization 
of the configurational distribution function piy,y) implies the 
following constraint for elements of the dual of the Lie algebra: 



6p(y, y)dy - 0. 



(22) 
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There is a further, more subtle constraint we need to impose 
on 5p{y,y). The configurational distribution function describes 
the probabihty for finding a shear deformation between any 
two layers, without considering any coupling between differ- 
ent pairs of layers. However, the total imposed shear strain has 
to be distributed over all the layers of the deformed crystal. We 
hence impose the further constraint that there is no freedom of 
varying the total strain. 



y Sp(y, y)dydy — 0. 



(23) 



We next need to define the constraints on the Lie algebra which 
has the physical interpretation of the space of observables (the 
pairing of the Lie algebra and its dual corresponds to the aver- 
aging of observables). In defining the functional derivative of 
some functional A of the configurational distribution function, 
the constraints (l22T l and (l23T l imply the following freedom: 

^ + f(y) + Cy, (24) 

5p(y,y) 6p(y,y) 

with an arbitrary function f(y) and a constant C, which play the 
role of Lagrange multipliers which still need too be fixed. In the 
spirit of reciprocal base vectors, where the dual of 1 is 1 and the 
dual of y is the differential operator d/dy, we use this freedom 
to satisfy the dual conditions 

6A 



/ 

md 



5p{y,y) 



dy^Q, 



6A 



-dydy - 0. 



(25) 



(26) 



dySpiy,y) 

Equations dZST l and ( |26] ) are the constraints on the elements of 
the underlying Lie algebra (of functional derivatives), whereas 
(l22l) and ( |23] ) are the corresponding constraints on the elements 
of the dual of the Lie algebra (of variations of the independent 
fields). The constraint (l26T l plays a crucial role in obtaining a 
meaningful evolution equation for the configurational distribu- 
tion function because it makes sure that the imposed macro- 
scopic shear deformation is accommodated by the crystal. Mo- 
tivated by the form of these constraints, we introduce two kinds 
of averages: (• • ■) for an average performed with the configura- 
tional distribution function p, and (■ • ■)sp for spatial averaging 
in the y direction. 

We can now evaluate the suitably normalized functional 
derivatives of energy and entropy (to be consistent with a three- 
dimensional description of the system one needs to introduce 
appropriate factors of area, Na^, into our formulation for unidi- 
rectional flows): 













^ Sp(y,y) ' \ 



vCy) 
1 



0(y)-<(l.>-(y-<r))((f)) 



,(27) 



and 















V 5p(y,y) ' 


V 





1/7^(3') 

-{Na\'kj,{\np{y,y)-{\np)] ) 



(28) 



where v{y) is the x-component of the velocity field and T{y) 
is the temperature field. Note that the last components in the 
functional derivatives (|27] ) and jTM satisfy the constraints dZSl) 
and ( l26b . 

From the general theory of Lie-Poisson reduction mentioned 
in the introduction, we obtain the following Poisson bracket: 



{A,B] = J J piy,y) ( 



d SA 



d 6B 
dy6p{y,y)j\dy6giy) 

d 6B \/ d 6A \ 
dy6p(y,y))\dy6g(y)l 



•dydy. 



(29) 



There are no convective contributions to this Poisson bracket 
because the system properties do not change in the direction 
of flow. The flow aff'ects only the configurational distribution 
function by changing the shear deformation y. The total en- 
tropy function is a degenerate function of this Poisson bracket 
because the configurational entropy density depends on y only 
through p. 

With the dissipative bracket, we express diffusion (or 
smoothing) in the space of shear deformations. 



[A,B] = 



1 



d 



p(r,y) 

6A d(^{y) 6A 



dy6p{y,y) 
d 5B 



dy 6e{y) 
d(b{y) 6B 



dy 6p(y, y) dy 6e(y) 



dydy. 



d(b\ 6A\ 
'djl 'S^l,, 
d<!)\ 6B\ 
'dyj '6^1, 

(30) 



where D is an inverse time scale for the difiiision of entire lay- 
ers, which we expect to be proportional to 1 /N. The energy is a 
degenerate functional of this positive-semidefinite, symmetric 
bracket. At this point, we have formulated all the thermody- 
namic building blocks of our toy model. We can now look at 
their implications. 

From the fundamental evolution equation (|6]l of the 
GENERIC framework iHIHIil], we obtain the following evo- 
lution equation for p: 



dp 



d_ 
dy 



(dv 

dy 



+ D 



1 £/<[> 

ksT dy 



D £/0' 



dy dy)'^^^^ 



According to GENERIC, in the momentum balance equation, 
there occurs the following elastic shear stress associated with 
the potential <1), 



\ ld<t> 



(32) 



From now on, we consider the GENERIC toy model for 
isothermal homogeneous shear flow, that is, the shear rate. 



r 



dv_ 

dy 



(33) 



and the temperature T are independent of the position y. As 
also the shear stress (l32l l becomes independent of y, the diff'u- 



sion equation dJTT i for the time-dependent configurational dis- 
tribution function p{-y) can be rewritten as 



dp 
5f 



d_ 

dy 



o-e D 



d I dp 
oy \ dy 



with the viscosity parameter 
fl3 ND' 



(34) 



(35) 



The anticipated dependence of the inverse time scale D on N 
mentioned after Eq. ( l30t has been built into this definition of 
fj; the parameter fj can hence be interpreted as the characteristic 
viscosity scale for single-atom processes and should hence be 
of the order of the viscosity of low-molecular- weight liquids. It 
is convenient to summarize the drift terms in equation ( l34b by 
introducing a potential with a constant external field. 



(l)(y) = (l)(y) - (fiy + a-e)Na^y. 



(36) 



The external field reflects the viscous and elastic stresses. 

From the diffusion equation ( l34l i for the configurational dis- 
tribution function, we obtain by averaging 



dt 



(37) 



By introducing constraints we have indeed forced the crystal 
to feel the full applied shear deformation. Given the integers 
Z as a discrete normal subgroup of the additive group of shear 
deformations leaving the potential energy invariant, Eq. (|5]l can 
be used to split the configurational shear deformation (y) into 
the elastic part 



1/2 1/2 

(r>e = (r>z = I ypz{y)dy = ^ I yp(z + y)dy. 



(38) 



-1/2 



1/2 



and the plastic part 
<r>p = <r) - <r>z = ^ zpz with p, = I p(z + y)dy. (39) 

zeZ ^ 



1/2 



-1/2 



The elastic part feels only the shear deformation with respect to 
a reference state; the plastic part is the shear associated with the 
reference states themselves. From the diff'usion equation ( l34l i. 
we obtain 



dt 



<r)p = f p 



r + — 

'7 



Pz. 



D 



dpz 
dy 



(40) 



y=I/2 



In the following, we consider only the diffusion equation (l34b 
rather than the complete set of thermodynamically admissible 
evolution equations. An important consequence of the general 
geometric ideas for formulating thermodynamically admissible 
models, however, is the fact that the argument y of the prob- 
ability density p should be recognized to represent the group 
of shear deformations R, and that the configurational energy is 
invariant under shear deformations from the discrete subgroup 



4.2. Solutions for s — 

The equilibrium solution (for y = and cTg - 0) of the diffu- 
sion equation ( l34l i for the potential is expected to be given 
by the Boltzmann factors. 



p'^'^iy) oc exp 



o.(y) 



(41) 



For a periodic potential, however, such a solution is not nor- 
malizable on the real axis R. We hence need a more careful 
discussion and we start with the potential Oo for e = 0, that is, 
with infinitely high energy barriers between wells. Then, one 
should note that for all integers z 



P(7)\ 



dpiy) 



r=z+ 1/2 



dy 



-0 for z e i 



(42) 



y=z+l/2 



SO that the points z+l/2 for z e Z cannot be passed or reached. 
We can hence look for separate solutions on each interval 



1 1 

2 2 



for z 6 Z 



(43) 



around any integer z- On any interval the solution (HTt can 
be normaUzed, 



/'"'(r) = 4 exp 



'i>o(y) 



1/2 

-1/2 



dy, (44) 



where the normalization constant Z is independent of z- The 
most general equilibrium solution on the entire real axis R is 
given by 



1>o(y) 



for ye/, with ^ p- - 1, (45) 



zeZ 



where the probability p, for a deformation in the interval /; does 
not change in time and is hence given by the initial conditions 
from which an equilibrium state is reached. For nonpatholog- 
ical initial conditions, also time-dependent solutions fulfilling 
Eq. (l42T i can be constructed independently on each interval /j 
and pieced together. 

We now apply a homogeneous shear flow that leads to a mod- 
ified force law (|36] | with a constant external force. For the par- 
ticular potential (B^ given in Eq. ( fT9] l, we obtain 



i>e(r) = '^siy) - if)y + o-^)Nay 



1 



2k^ (cos ny)^ + s 



- ( fir + o"e)r 



(46) 



For e = 0, we can repeat the previous construction of steady 
state solutions fulfilling the consistency conditions (l42l) for any 
finite applied shear rate y. The solution (l44l) on the interval 
is replaced by 



1 / i>()(r) 



exp 
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<E'o(y) 



dy, (47) 



where the normalization constant Z, now depends on z- The dis- 
tribution on each interval L is deformed by the flow, but there 
are no transitions between different potential wells. According 
to Eq. (|39] l. we hence have - 0. Under steady-state condi- 
tions, we have also - 0. 

4.3. Solutions for small e 

If we proceed from e = to small values of e, the barriers at 
z + Ijl become large but finite so that there arises the possibil- 
ity of transitions between neighboring intervals. The transition 
rates are small so that the probabilities p~ change slowly. The 
evolution in each interval is much faster so that a steady state 
is reached in each interval before there is a substantial proba- 
bility for transitions. To the next level of approximation, the 
solutions of the diffusion equation (l34l i are hence still given by 
Eq. (l45T l. but the probabilities p- are to be obtained from a set 
of evolution equations in terms of the (small) transition rates. 

Let us consider transitions between the intervals L and in 
the absence of flow (for j - Q and cr^ - 0), when the probability 
p. is different from p,+i. In view of this diff'erence, the almost 
stationary solutions around the peak at z H- 1 /2 cannot be of the 
flux-free symmetric Boltzmann form (HTI) . We need a skewed 
stationary solution that can interpolate between two different 
levels of probability. The solution to this problem is given by 
the following generalization of the Boltzmann distribution in 
Eq. ( |4TI ). which is fundamental for the subsequent analysis: 



p'^ij) = exp i - 
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(48) 



where c and c' are suitable constants fixing the normalization 
and the flux. Note that c depends on the choice of the lower 
limit of integration, yo, whereas c' does not; for small c', how- 
ever, this dependence is weak. For the solution ( |48] ), the proba- 
bility flux is given by the constant 



-D\ -p"-^ + -^\^ c'D, 

k-bT ay ay 



(49) 



so that we obtain a stationary solution of the diffusion equation 
( l34b with a constant probability flux through the system. By 
using Eq. ( |48] | around z+ 1/2, we obtain for small differences 
between p, and p,^i 
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where the integral 
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when evaluated with the saddle point approximation, is inde- 
pendent of the precise value of the small parameter 6 (with 
e «; 5 <K 1) because it is dominated by the high potential bar- 
rier around y - 1/2. Equation (l50t determines the small flux 
or flow rate parameter c'. With this value of c' and c = Pz/Z, 



we obtain the flux of probability from L to L+i. Taking into 
account the change of p^ also from the flux of probability from 
to the equations governing the slow time evolution of p^ 
can then be written as 



^ = riPz+i - 2pz + Pz-\), 
dt 

where the rate parameter r is given by 
_ D 

zi' 



(52) 
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In the presence of shear flow, matching solutions on neigh- 
boring intervals can be constructed in the same way. A compre- 
hensive solution summarizing all the features to the discussed 
level of precision can be written as 



p(z + y) 
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■ exp 
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exp \ [ dy' for y e [-1/2, 0] 



on the interval I- (note that the solution on I- is conveniently 
shifted to the interval /q). The definition of the integral / of 
Eq. ( BTT l remains unchanged and the parameter w is given by 



(55) 



The integral terms in Eq. ( l54l i are negligible except in narrow 
boundary layers. The values of the solution ( |54] | on the bound- 
aries are given by 
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so that the continuous matching of the solution on neighboring 
intervals becomes obvious. Moreover, we obtain 
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The modified time-evolution equations for p- in the presence 
of flow become 



dPz 
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From Eqs. (|39] l and ( |59] l, we obtain 
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(60) 



which is consistent with Eqs. ( |40] | and ( ISTl i. Equation ( |60] | re- 
lates the plastic strain rate to the externally applied stress 
(note that fjj is expected to be small compared to o-g). It is 
closely related to Eyring's viscosity formula fiB. The "vol- 
ume of a flow unit" or "viscosity volume" is given by Na^, 
where such a large volume is a hallmark of highly viscous, plas- 
tic flow. Accordingly, the rate parameter ro is extremely small. 
Equation ( l60l l is the main result of the asymptotic solution of 
our toy model on the full universal covering group. Equation 
dSTb shows that the configurational probability density on the 
quotient group occurs naturally. The evolution equation ( |59] ) 
implies the following H-theorem-like identity for transitions be- 
tween the different intervals. 
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For a weak dependence of p, on z, Eq. 
as a partial differential equation. 
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< 0. (61) 
can be rewritten 



(62) 



We once more realize that, on top of slow diff'usion, there is a 
small systematic flux of probability with the drift velocity j-^ 
given in Eq. (|60] |. 

5. Toy model on quotient group 

So far, we have developed a kinetic theory on the full, sim- 
ply connected group of shear deformations. Equation ( |59] | de- 
scribes the irreversible transitions between reference states. In 
particular, it implies the plastic deformation rate (l60l l. We now 
address the construction of a complete kinetic theory for the 
relative deformations, that is, on the quotient group Iq - R/Z. 
We look for an approach that reproduces the results of the pre- 
ceding section, but works directly and exclusively with a con- 
figurational distribution function on the interval Iq. 

5.1. Key equations 

The general definition (|5]i of a probability density on Iq = 
R/Z for the group of shear deformations reads 



Pzir) = Z /'(z + r) for r e k. 



(63) 



as already introduced in Eq. ( l38l l. Because of the periodicity 
of the potential and the linearity of the diffusion equation 
(O in p for given cTg, from any solution p on the entire real 
axis R, one obtains a solution pz of the diffusion equation on 
Iq. If the original solution is smooth, this construction ensures 
that we obtain the same smoothness for the periodic boundary 
conditions because one has identical values at y = -1/2 and 
J - 1/2 in Eq. ( l63T l, as can be seen by shifting the summation 
index z- A kinetic theory on /q must hence be based on smooth 
periodic solutions of the diffusion equation, so that we actually 



deal with a circle. Note that also the elastic stress o-f. given in 
Eq. ( l32l l can be evaluated as an average performed with pz, 



1 /c/ct 



A^fl^ \ dy 



(64) 



because the averaged force is a periodic function of the shear 
deformation. 

In summary, we need to solve the following diffusion equa- 
tion with periodic boundary conditions on Iq: 



dpz 
dt 
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From the diffusion equation (1651 ), we obtain by averaging 



:T-<r>z = r- 

dt 



7+ — \px-D— 
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= r-Tp- 



(66) 



y=i/2 



We thus naturally recover the splitting of the total shear rate 
into an elastic contribution and the plastic contribution ( l40l i. 
The model with a configurational distribution function p on the 
real line hence justifies the following assumptions of the kinetic 
theory with pz on the interval Iq: (i) the use of smooth periodic 
boundary conditions for pz. (ii) the shear stress CTg given by 
Eq. ( l64l i. and (iii) the splitting of the shear rate into elastic and 
plastic contributions, where the latter is given in Eq. ( l40b . All 
these assumptions have been made in the original formulation 
of a toy model of crystal viscoplasticity ifioll . Formally, the peri- 
odic theory works perfectly well for time-dependent problems, 
even when the energy barrier between wells is not assumed to 
be very large. 



5.2. Steady state solution 

For steady state situations, the solution of the diffusion equa- 
tion ( 1651 ) is of the general form given in Eq. ( l48b . The constants 
have to be chosen such that we obtain a periodic normalized 
solution. In particular, we can then determine the plastic shear 
rate from Eq. ( |66] | tlOll : 



Jexp{- 
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-1/2 



dy' 



dy 



Dlw 



(67) 



For small e, this result is consistent with Eqs. (fSSl l and (|60] |. A 
Kramers-type formula fl3*l for the rate of transitions between 
reference states thus follows immediately from the periodicity 
requirement. In some analytical developments, but certainly 
also in numerical calculations, it is more convenient to obtain a 
solution of the diffusion equation on a compact interval with pe- 
riodic boundary conditions rather than a drifting and diffusively 
broadening solution on the real axis. 
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5.3. Thermodynamic formulation 

At this point we know that a complete and elegant kinetic 
theory can be constructed on the quotient group of relative de- 
formations with respect to a slowly evolving reference state. 
However, several further questions remain to be addressed: Can 
one find thermodynamic building blocks directly for the quo- 
tient group model? If so, how are these building blocks related 
to those of the model on the cover group? Are there any ex- 
tra rules or recipes that need to be added to the thermodynamic 
framework? 

An obvious additional recipe is the periodicity requirement 
for p-z- This periodicity cannot be derived directly from the 
thermodynamic consistency on the cover group; rather, peri- 
odicity is a consequence of superposing all smoothly patched 
solutions on the different equivalent intervals. In that sense, the 
periodicity requirement is a fundamentally new ingredient of 
the thermodynamic approach resulting from the act of cover- 
ing. For generalizing the periodicity requirement to more com- 
plicated situations, it is important to remember that there exists 
a well-developed theory of harmonic analysis on locally com- 
pact groups (with particularly powerful results for commutative 
or compact groups). 

We now turn to the thermodynamic building blocks. The 
conformational contribution to the energy in Eq. ( fTSl l per unit 
channel width, (l/fl)(((I)))sp, can be rewritten as (l/fl)((0)z)sp 
because the potential energy is invariant under shear transfor- 
mations from the discrete group Z. The discussion of entropy 
can be based on the following identity: 

oo 

— CO 

+ /jP^i — 1" : T^^- ^^^^ 

^ Jl Pz PzPz(7-z) 

The first term on the right-hand side of Eq. ( l68T l represents the 
entropy associated with the periodic probability density pz on 
the interval /q, whereas the second term results from the distri- 
bution over the equivalent intervals The last term vanishes if 
the normahzed probability on each interval /j, that is p{y)lpz, 
is equal to the properly shifted distribution pz- Any deviation 
of piy)/p: from the shifted pz leads to a reduction of entropy. 
This loss of entropy cannot be estimated with the knowledge 
of Pz and alone. Nevertheless, the dynamics on the inter- 
val Iq can be described rigorously. Therefore, the entropy term 
resulting from a lack of proportionality of the distributions on 
the different intervals /, can only affect the transitions between 
the different intervals. The transition regimes between different 
intervals are not self-similar and cannot be obtained from the 
theory on the quotient group. 

After discussing the conformational contributions to the en- 
ergy and entropy in terms of pz, we need to calculate their func- 
tional derivatives. In particular, we need to analyze the signif- 
icance of the constraint ( l23T l on the average shear deformation 
and the reciprocal condition (l26b in defining functional deriva- 
tives for the quotient group model. The constraint condition 



( 1261 ) was recognized to be crucial to obtain the proper diffu- 
sion equation; we hence need to take it over into the model on 
the quotient group. In calculating dEjSp, this constraint leads 
to a jump of SE/dp rather than periodic boundary conditions. 
Such a jump is actually crucial to obtain periodic steady-state 
solutions pz in the presence of a flux. Some overall effect of 
the jumps in the probability levels p, of the smoothly patched 
solutions on the different intervals /, due to a flux in a driven 
system, which is lost in defining pz by superposition, is reintro- 
duced through the functional derivatives with jumps. Note that 
the functional derivatives are modified such that there occurs a 
jump between y = -1/2 and y - +1/2, but the periodicity of 
all derivatives with respect to y remains unaffected. The proper 
definition of functional derivatives has to be seen together with 
the new periodicity requirement. In the original direct construc- 
tion of a kinetic theory on the interval Iq [10] , the counterpart of 
the constrained functional derivative is achieved through a for- 
mal driving parameter which can be related to the macroscopic 
strain rate and shear stress by imposing the constraint that the 
applied strain rate is the sum of the elastic and plastic strain 
rates. 

We finally need to discuss the Poisson and dissipative brack- 
ets. The Poisson bracket results from a Lie algebra and hence 
coincides for the group Iq and its cover R. Also the dissipa- 
tive bracket in Eq. dSOl l involves only local quantities and an 
average of a periodic function, so that it may be used directly 
on the quotient group. When all the thermodynamic building 
blocks are used in the fundamental equation we recover the 
equations of Sect. 15.11 Periodicity supplemented by properly 
defined functional derivatives are the new ingredients into the 
thermodynamic framework. 

6. Summary and discussion 

The present paper can be considered from two different per- 
spectives. From a mathematical perspective, it analyses the 
freedom in associating Lie groups with the Lie algebras that 
provide the structure required to introduce reversible time evo- 
lution. From a thermodynamic perspective, this paper develops 
the tools for models with multiple equivalent reference states 
evolving slowly in time. 

All continuous groups obtained as quotients of a universal 
covering group and a discrete normal subgroup of the cover- 
ing group possess the same Lie algebra. The discrete normal 
subgroup, which may be finite or infinite, may be interpreted 
as a set of equivalent reference states, and the quotient group 
describes deformation effects with respect to a reference state. 
For a toy model of crystal viscoplasticity, we have shown in de- 
tail how modeling on the cover and quotient groups works, and 
how they are related. On quotient groups, the new ingredients to 
thermodynamic modeling are suitable periodicity requirements 
supplemented by properly defined functional derivatives. Im- 
portant mathematical ingredients in addition to the correspon- 
dence between Lie groups and Lie algebras are Haar measures 
and harmonic analysis on locally compact groups. 

The ideas of this paper naturally lead to a classification of 
rheological models, depending on whether or not we make use 
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Figure 1: Attributes of theological models. 



of reference states (see Fig.[T]). From an abstract point of view, 
reference states occur when the Lie group providing the struc- 
ture required to set up thermodynamics is not simply connected 
so that it can be written as the quotient of the universal cover- 
ing group and a discrete group of reference states. In partic- 
ular, we obtain a clear distinction between viscoplasticity and 
viscoelasticity, which often is not very transparent in the liter- 
ature, in terms of reference states. In the absence of reference 
states we obtain the models of viscoelasticity which have been 
applied successfully to many complex fluids. The possibility of 
transitions between equivalent reference states is the hallmark 
of viscoplasticity. 

Viscoplasticity is often defined as rate-dependent plasticity, 
where plasticity generally means the ability to change or de- 
form permanently (in contrast to elasticity, which refers to the 
ability to change temporarily and revert back to the original 
form). These permanent changes occur by transitions from ref- 
erence states to equivalent ones. Within the domain of a par- 
ticular reference state, deformations are reversible and hence 
elastic. Only when a sufficiently high stress level is reached, 
transitions between equivalent reference states occur with a no- 
ticeable rate (in principle, transitions are possible even in the 
absence of stresses, but they are extraordinarily rare). This 
stress level is known as the yield point or, for anisotropic ma- 
terials, as the yield surface. As yield is a matter of "notice- 
ability," phenomenological definitions of viscoplasticity based 
on this concept may not be ideal. The deformation behavior 
described by long-living reference states is often referred to a 
elastic-viscoplastic, whereas transitions at a noticeable rate lead 
to a plastic flow associated with a very large viscosity. 

The deformation with respect to evolving reference states en- 
ters many models of viscoplasticity in which the total defor- 
mation is split into a plastic and an elastic part, where only 
the elastic part is used as a structural variable (see Leonov 
131, Hiitter and Tervoort fll'] and references therein). Such a 
splitting corresponds directly to the factorization of a covering 
group. Usually, the relative or elastic deformation is treated by 
explicit variables, whereas the transition rate between reference 
states merely contributes to the rate of change of the relative 
deformation variable. The relative deformation measures used 
in previous phenomenological models of viscoplasticity corre- 



spond to moments of the configurational distribution function 
in the kinetic theory of the present paper, just as the conforma- 
tion tensors of the Maxwell and Oldroyd models of viscoelas- 
ticity correspond to second moments of the connector vector of 
a Hookean dumbbell model in polymer kinetic theory 121]. Even 
if a kinetic theory on the quotient group does not treat reference 
states explicitly, the suggested analysis of the group theoretical 
background clearly shows the underlying existence of reference 
states and hence the viscoplastic nature of the models. Group 
theory thus helps to identify viscoplastic models, and the phe- 
nomena described by such models in a meaningful way are then 
recognized as viscoplastic. 

Kinetic theories of viscoplasticity, such as the toy model em- 
ployed for the purpose of illustration, are a promising starting 
point for deriving existing phenomenological models and for 
generalizing them, both for crystalline and for amorphous ma- 
terials. They offer a link between detailed atomistic considera- 
tions and phenomenological models 127, 22, "s*]. For 
deriving kinetic theories by statistical mechanics 119, 
Eoll . it is important to realize the occurrence of discrete sets of 
reference states and of transition rates. The statistical mechan- 
ics of models with reference states should hence be expected to 
be based on Kramers-type formulas lll3r for transition rates. 
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